Non stationary wave turbulence in an elastic plate 
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We report experimental results on the decay of wave turbulence in an elastic plate obtained by 
stopping the forcing from a stationary turbulent state. In the stationary case, the forcing is seen 
to induce some anisotropy and a spectrum in disagreement with the weak turbulence theory. After 
stopping the forcing, almost perfect isotropy is restored. The decay of energy is self similar and the 
observed decaying spectrum is in better agreement with the prediction of the weak turbulence theory. 
The dissipative part of the spectrum is partially consistent with the theoretical prediction based on 
previous work by Kolmakov. This suggest that the non agreement with the weak turbulence theory 
is mostly due to a spurious effect of the forcing related to the finite size of the system. 
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PACS numbers: 46.40.-f,62.30.+d,05.45.-a 

Wave turbulence is observed in systems involving a 
large number of waves coupled by nonlinear effects. Like 
hydrodynamical turbulence, it is characterized by wide 
ranges of excited length scales and time scales. As wave 
turbulence can be weakly nonlinear, a statistical theory 
(Weak Turbulence Theory, WTT) could be analytically 
derived and applied to numerous systems such as optics, 
oceanic/atmospheric waves or plasmas (see [l|, [^l for re- 
views) . The predicted phenomenology of the forced case 
is very similar to that of hydrodynamical turbulence: the 
energy injected at large scales where dissipation is negligi- 
ble is transferred conservatively to small scales by nonlin- 
earity through the Kolmogorov-Zakharov (KZ) cascade. 
The cascade operates until the energy reaches scales that 
are small enough so that dissipation dominates over non 
linear transfers. In this respect wave turbulence can ap- 
pear as an intermediate step in complexity towards the 
understanding of hydrodynamical turbulence. 

Assuming that the waves composing the motion ex- 
change energy at a slow rate compared to their frequency 
(the non linear interactions are weak) and that the sys- 
tem is infinite, the WTT yields analytical predictions 
for spectra and other other statistical quantities. Let 
flk stand for the canonical variables that diagonalize the 
linear part of the Hamiltonian of the motion. The wave 
action (or occupation number at wavevector k) is defined 
as TT-k — {aiti{t)a^(t)), where () is a statistical average on 
realizations and * stands for complex conjugaison (see [3[ 
for details in the case of flexion waves on a thin elastic 
plate). In the framework of the WTT, the time evolution 
of rik can be described by 
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The forcing Fk is usually acting on large scales and the 
dissipation ^\ji\^ is efficient at small scales. In the in- 
termediate inertial range of scales, energy is transferred 
conservatively in Fourier space. Col{n\^) is the collision 
integral whose expression was derived in many physical 
systems [l|-l3| . It appears as the divergence of an energy 
flux in Fourier space and formally it is reminiscent of the 



kinetic theory of gases as the energy is transferred non- 
linear ly by "collisions" of resonant wavetrains [l|, l2| . The 
Kolmogorov-Zakharov (KZ) spectrum is a stationary so- 
lution of the forced case whose analytical expression can 
be derived in the intermediate range of scales. 

Here we consider the case of flexion waves in an elastic 
plate. In the experiment we measure the velocity spec- 
trum £'(k) which is related to rik by E{k) oc /c^rikQ]. 
The theoretical prediction for the velocity spectrum is: 
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where C is a dimensional constant that can in principle 
be calculated analytically, P is the average input power 
and fee is a cutoff frequency that has to be introduced for 
the sake of dimensionality . The logarithmic dependency 
has to be introduced to ensure a flnite flux of energy be- 
cause of the degeneracy with the equipartition spectrum 
-E(k) — constant which has a zero flux of energy [3]. The 
relation between the cutoff wavenumber fcc and the small 
(or large) scale dissipation is largely an open question 
as no dissipation is considered in the theory. Although 
Diiring et al. claim to have observed the KZ spectrum in 
numerical simulations (by choosing kc as the dissipation 
scale of the numerical simulation) so far the KZ spectrum 
remains elusive in experiments as the theoretical scaling 
in P and in k (or w) is not observed [J, |5|. More gener- 
ally, the applicability of the Weak Turbulence Theory to 
real systems remains an open question and only few ex- 
perimental systems are available for a quantitative com- 
parison with the theoretical predictions (see the recent 
review [l[ for a detailed discussion). In this respect, the 
case of flexion waves in a thin plate is promising because 
the use of a high speed profllometry technique enables 
advanced measurements of the full space-time structure 
of the turbulent fleld 5-7] . Although not in agreement 
with the WTT, the flexion wave turbulence has been ob- 
served to be indeed weakly non linear, with a persistence 
of waves and a non linear dispersion relation weakly mod- 
ified from the linear dispersion relation. These observa- 
tions are consistent with the requirement and the pre- 
dictions of the WTT. It has been suggested that either 
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dissipation or finite size effects could be responsible for 
such discrepancy [j-|7|. Here we study a non stationary 
regime of such turbulence obtained by stopping the large 
scale forcing after a stationary regime is reached. By 
observing the decline of turbulence, we suggest that the 
forcing itself is responsible for the disagreement between 
the measured spectrum and the theoretical one through 
finite size effects. 

Non stationary weak turbulence was studied theoret- 
ically or numerically by several authors [2, IsMld] but 
not in conditions that are relevant to laboratory stud- 
ies. Usually dissipation is not considered and energy is 
initially localized in a finite interval of (large) scales so 
that the propagation of an energy front in Fourier space 
to large wavenumbers was observed [9|, llO| . The experi- 
mental situation of decaying wave turbulence where dissi- 
pation is present and starting from stationary turbulence 
is different and was studied experimentally for capillary 
wave turbulence ll[ but with limited measurements due 
to the difficult experimental conditions. A self similar 
decay of the spectra was predicted in this case 
has not been observed yet. 
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FIG. 1: Schematics of the experimental setup (similar to [a|). 
A 2 X 1 m^ plate, 0.4 mm thick, is vibrated by an electromag- 
netic shaker at 30 Hz. A single value of the forcing amplitude 
is studied here. The deformation of the plate is measured 
by a high speed profilometry technique [l3l | over a surface 
1.25 X 0.94 m^. Linear waves follow a dispersion relation 
00 = ck^ with c close to 0.64 m^.s~^ so that the 30 Hz forc- 
ing corresponds to a wavelength A ~ 0.37 m. Movies of the 
deformations were recorded at 6000 frames/s. 

The forcing is stopped after the plate reached a sta- 
tionary turbulent state. Averages over realizations are 
performed by repeating the experiment 90 times (see 
fig. [U for a sketch of the experimental setup) . Movies of 
the deformation of the plate are obtained by high speed 
profilometry [6|. Fourier analysis of the data is similar 
to [5[. At a given k, two counter propagating waves are 
present with frequencies of opposite signs so that the spa- 
tial Fourier transform of the velocity field can be written 
generally as u(k, t) = u+(k, t) + v~ (k, t) where w"*" (resp. 
v~) contains only positive (resp. negative) frequencies. 
We study only Fourier components of negative frequen- 
cies (i.e. v^) by applying a discrete Hilbert transform 
in time (without loosing any information as v^ (k, t) — 
v~*{—\i.,t)). In this way, we can separate waves at a 
given k traveling in the two opposite directions so that 
to study isotropy issues. We chose to perform the statisti- 
cal analysis on v~ so that for instance waves with fc^; > 




are propagating in the direction a; > which makes the 
discussion slightly easier. Here we define the 2D time 
dependent Fourier spectrum as -E(k,t) = (|u^(k, t)|^). 
Note that this spectrum is no longer symmetrical under 
the change k — > — k. E{k = l|kl|,t) is the angle averaged 
spectrum E{k,t) — J^^ E{'k,t)kd9 where 9 is the polar 
angle of k. To get some insight on anisotropy issues, 
the interval of integration can be restricted to four quad- 
rants of angles in [-7r/4,7r/4], [7r/4,37r/4], [37r/4, 57r/4] 
and [57r/4, 77r/4]. These intervals correspond to propa- 
gation along directions toward x>0,y>0,a;<0, 
y < respectively (noted a;+, y+, x" , y^ resp. in the 
following) . 
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FIG. 2; (a) Decay of the velocity spectrum E{k — 2tt/X, t) 
at various times. Top solid curve (S): stationary spectrum. 
Lower solid curves, from top to bottom: 104 ms after stop- 
ping the forcing and every 83 ms subsequently (spectra are 
averaged over 42 ms). Thin dashed line at the top of the fig- 
ure: shape of the theoretical spectrum (_B(fc) oc fcln^'^(fcc/fc)). 
Thick dashed line: fit of the model ([3]) for the dissipative part 
of the spectrum at each time. Inseri: rescaling of the decay- 
ing spectra using the fitting parameter of Q at each time as 
a function of k/kb{t) (covering both the inertial and the dissi- 
pative range) for times between 0.1 and 0.83 s. Upper dashed 
line: power law fc"'^, lower dashed line: theoretical spectrum, 
(b) fitting parameter B /k1{t) as a function of time. Dashed 
line: linear trend, (c) fitting parameter A as a function of 
time. 

An example of the decay of the velocity spectrum is 



shown in fig. [UJa). The top sohd curve shows the sta- 
tionary spectrum observed before stopping the forcing 
shaker. A maximum is observed at the forcing wavenum- 
ber and a short power law inertial regime is observed with 
an exponent about —0.2 corresponding to the previously 
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scaling of the single point frequency spec- 



trum (change of variable w -O- fc is obtain through the 
dispersion relation cj oc k"^) [J, LZ|- The theoretical spec- 
trum (121) is shown with a thick dashed line and is clearly 
different from the experimental one both in the fact that 
the latter is increasing in the inertial range and that the 
cutoff is not following a logarithmic decay. After stop- 
ping the forcing and after a short transient during which 
the forcing peak disappears, the spectra are decaying in 
a self similar way until the inertial range is suppressed. 
A peak at low k is observed which corresponds to the 
energy of the lowest frequency modes of the plate which 
are not resolved by our measurement (the field of view 
being too small). The cutoff wavenumber decays with 
time. The exponent of the inertial range scaling is ob- 
served to be close to 0.3 i.e. closer to the theoretical 
exponent (equal to one if one discards the log correction 
in E{k) ex fcln^'"^(fcc/fc)) than for the forced case (—0.2). 
The remaining disagreement may be attributed to a too 
narrow interval of length scale available for the energy 
cascade or to the logarithmic corrections. 
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FIG. 3: Decay of energy. Main figure: evolution of the angle 
averaged power spectrum E{k, t) for 1/A = fc/27r = 5, 10, 15, 
20, 25 m"^ from top to bottom (semilog scale, the curves have 
been shifted vertically for clarity). Dashed lines: eye guides 
for the long time exponential decay of the energy. Insert: evo- 
lution of exponential decay times as a function of k. Dashed 
line: Lorentzian fit (0.73 + 0.025(fe/27r)^)"\ 

The time evolution of the energy at a given k is shown 
in figure [3] for a few values of k. The decay occurs in 
two steps: the first one resembles a stretched exponen- 
tial shape, the second one is exponentially decaying. The 
change of behavior in the decay is related to a change in 
dynamics from the non linear cascade to the dissipative 
decay. The characteristic times of the exponential decay 
are shown in the insert of fig. [31 For k above the forcing 
scale, they tend to follow a Lorentzian variation with k 



suggesting the coexistence of two processes for dissipat- 
ing energy at large or small scale (the increase at the 
upper values of k is due to a very poor signal over noise 
ratio). The values of the dissipative times (from 1 s to 
1/10 s) show also a clear scale separation from the pe- 
riod of the wave: the ratio equals 25 at the forcing scale 
up to over 100 at k/2'K = 25 m~^. The exponential de- 
cay is consistent with the dissipative term —jkn^ with 
7fc = 0.73 + 0.025(fc/27r)2 s^^. The observed time scale 
separation opens the possibility of an energy cascade as 
predicted by WTT. 

A theoretical prediction for the shape of the decaying 
spectrum can be derived following the theoretical work 
of Kolmakov [l^: one looks for a self similar decaying 
solution of ^ expressed as nk{t) = Akt{t)°'g{£,) where 
kbit) is a time dependent dissipative cutoff wavenumber 
and £, — k/kh{t). We assume that 7k ex fc^ as suggested 
by the data. In order to express ([1]) as a function of 
the single parameter ^ for a self similar decay of the 
spectrum, the scaling properties in k of the collision in- 
tegral (as given in [3|) and of the dissipation time im- 
pose that a = in the case of the plate and the cutoff 
wavenumber must verify kh{t)kb{t)~^ = constant. One 
gets kl{t) = fc^(0)/(l + t/r) where r is a constant char- 
acteristic time. In the dissipative range k ^ ki,(t), the 
collision integral in ([T|) is negligible compared to the dis- 
sipative term so that ^ yields g'{S,) = —2BS_g{^) (where 
the prime stands for derivation and _B is a positive con- 
stant) so that g{^) ex exp(— B^-^). The overall dissipative 
region of the spectrum is thus expected to follow 



Edik.t) = Ak^e^p{-Bk^/kl{t)) 
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where A is a positive constant depending on the spectrum 
at initial time. The dissipative cutoff is then expected to 
follow a Gaussian decay with l/k'^{t) being affine with 
time. The decay ^ has been fitted to the dissipative 
region in fig. [2] that provides A and B/k'^(t) (fig. [2] (b) 
and (c)). The agreement is fairly good and the cutoff 
wavenumber follows the theoretical prediction (l/fc^(t) 
being affine with t). Nevertheless the prefactor A of the 
exponential decay is not constant as predicted above but 
it is strongly decaying with time. Thus, the shape of the 
spectrum is preserved but a genuine self-similarity is not 
observed. It is interesting to notice that the prediction 
of kb{t) and the dissipative behavior of the spectrum ac- 
tually depend on the scaling properties of the collision 
integral C{n\^) and that they are at least partially com- 
patible with the data. Note that the cutoff observed in 
the stationary forced regime is also close to such a Gaus- 
sian decay. In the inertial range {k < kb{t)), the collision 
integral dominates and the spectrum is expected to follow 
the KZ spectrum |12| . 

The issue of isotropy is addressed by separating the 
fully angle averaged spectrum E{k, t) into the four quad- 
rants defined above (as shown in fig. [5]). In the forced 
case, waves in the quadrant x'^ dominates strongly. Due 
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FIG. 4: Check of isotropy of the spectrum, (a) stationary 
forced regime, (b) decaying case after 0.18 s of decay. The 
upper thick curve is the fully angle averaged spectrum. The 4 
lower curves correspond to an angle average over the 4 quad- 
rants (see text). Curve 1 corresponds to propagation direc- 
tions x^ , curve 3 to x~ and 2 and 4 to both directions along 
y (almost superimposed). 



to the position of the measurement region, the spherical 
waves directly generated by the shaker are only partially 
visible and lie mostly in this x+ quadrant. This explains 
that the energy in this quadrant is dominating over the 
other quadrants at the forcing wave number and most 
of the inertial range. Waves propagating in the other 
quadrants have been generated by nonlinearities and/or 
rebounds on the boundaries of the plate. The forcing is 
only weakly visible on the curves 2 and 4 (y+ and y" 
quadrants ) and not visible in the x~ quadrant . In con- 
trast to the forced case, the decaying spectrum is almost 
perfectly isotropic. Note that the peak at the lowest k 
corresponds to energy in the lowest frequency modes with 
wavelengths larger than the field of view of our technique. 
Some energy is stored in these modes although their fre- 
quency is lower than that of the forcing. Their energy 
is also decaying but may act as an energy reservoir that 
feeds the cascade during the decay. It should be noted 
that the spectrum of the x~ quadrant is not much differ- 
ent in the forced and the freely decaying regimes. 

We interpret our observations as follows: In the weakly 
forced regime, non linear effects are due to cumulative 
effects over several periods of the waves. Due to the 
propagating nature of the waves, it translates into some 
propagation length required for the generation of other 
frequencies by weak non linearities. If this distance is 
not short compared to the size of the plate, the transi- 
tory regime dominates: this is the case in the x~^ quad- 
rant. By contrast the waves in the x~ quadrant have 
propagated further (at least one rebound) and this seems 
enough for non linear effects to have transferred energy 
from the forcing to other waves. This explains why the 
spectrum of the waves propagating in this direction are 
close to that of freely decaying turbulence. We expect 
that a bigger system and a measurement region further 



from the forcing would yield a better agreement with the 
theory. Thus our observation can be assigned to finite 
size effect, the measurement region being too close to the 
forcing point. It is very different from another kind of 
finite size effect related to the quantization of the plate 
modes due to its size. The latter effect is expected to 
influence the energy cascade by limiting the number of 
wave resonances available to transfer energy 14 1 . 

The exponent of the isotropic decaying spectrum 
is closer to the WTT prediction of the Kolmogorov- 
Zakharov spectrum but not equal. The remaining dis- 
crepancy can be attributed either to a too short inertial 
range and/or to the logarithmic correction predicted by 
WTT. The self similar decay implies that the logarith- 
mic cutoff kc should evolve in time following the decay of 
kb{t) and thus display a dissipative scaling (as implicitely 
assumed in the numerical simulations of [3||). The par- 
tial agreement of the observed and predicted dissipative 
cutoff provides some support for the validity of the ex- 
pression of the collision integral derived by the WTT. 
The fact that the prefactor A is not constant in time is 
most likely due to the fact that the damping coefficient 
7fc is not purely quadratic in k. The constant part of 
7fe may accelerate the decay of energy as dissipation oc- 
curs at all wavenumbers. This may be responsible for the 
observed decay of A. 
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